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Abstract 

We formulate a low energy effective theory describing phases of matter that are both solid and 
superfluid. These systems simultaneously break translational symmetry and the phase symmetry 
associated with particle number. The symmetries restrict the combinations of terms that can 
appear in the effective action and the lowest order terms featuring equal number of derivatives 
and Goldstone fields are completely specified by the thermodynamic free energy, or equivalently 
by the long-wavelength limit of static correlation functions in the ground state. We show that the 
underlying interaction between particles that constitute the lattice and the superfluid gives rise to 
entrainment, and mixing between the Goldstone modes. As a concrete example we discuss the low 
energy theory for the inner crust of a neutron star, where a lattice of ionized nuclei coexists with a 
neutron superfluid. 



1 



1 Introduction 



The low energy dynamics of strongly interacting solids and superfiuids can be systematically studied 
through an effective theory formulation in terms of weakly interacting phonons - the collective degrees 
of freedom in these systems. In the familiar case of solids, one longitudinal phonon and two transverse 
phonons arise as Goldstone modes due to the breaking of translation symmetry. In the case of a 
superfluid, one mode called the supcrfhiid phonon arises due to the breaking of the global U(l) symmetry 
associated with phase rotations of a field operator 0. In special cases the ground state of the system 
can spontaneously break both these symmetries. A particularly simple but non-trivial realization is a 
solid immersed in a superfluid with strong interactions between the particles that form the solid and the 
superfluid respectively. It is likely that a substantial region in the crust of a neutron star is occupied by 
such a phase [1] and its presence may affect neutron star phenomenology. From general considerations 
we can argue that the inner crust of neutron stars features a lattice of neutron rich nuclei in a bath 
of unbound superfluid neutrons. The lattice sites can be viewed as clusters of protons, with a fraction 
of neutrons "entrained" on the clusters [5J [3] . Other intrinsically more complex phases where a single 
component exhibits both superfluid and solid characteristics have also been proposed. They include the 
supersolid phase of 4 He [4] and the Larkin Ovchinnikov Fulde Ferrell (LOFF) phases [5j [6] in polarized 
fermion superfiuids. Although these systems can in principle be realized terrestrially, they have proven 
to be challenging to explore in experiments [7]. Nonetheless in all these cases the low energy dynamics 
is described by an effective theory of four Goldstone modes [8]. The associated fields for the lattice 
phonons are £° =1 " 3 (r, t) and are related to space-time dependent deformations of the lattice. Similarly, 
the field associated with the superfluid mode <fi(r, t) is related to the space-time dependent phase of the 
condensate. Because of interactions, such as those between the neutrons and the protons in the neutron 
star crust, one can not in general treat the two sectors separately and a unified treatment is required. 
It is the aim of this paper to provide such a framework. 

The low energy theory is described in terms of the fields and £ a . The symmetries associated 
with translation and number conservation require that the low energy theory be invariant under the 
transformation £ a=1 - 3 (r,i) — > £ a=1 " 3 (r,i) + a a=1 - 3 and (f>{r,t) — > (p(r,t) + 9 where a° =1 " 3 and are 
constant shifts. This naturally implies that the low energy lagrangian can contain only spatial and 
temporal gradients of these fields. Further, by requiring cubic symmetry for the crystalline state, the 
quadratic part of the effective lagrangian is given by, 

£ = f (5 O 0) 2 - ^W) 2 + fordo? - \»(e b e b ) - fanm 6 ) 

o=1..3 

where higher order terms involve higher powers of the gradients of these fields, and £ ah = (<9 a £ b + <9&£ a ) — 
|9 c £ c <5 ah . In the uncoupled case, the low energy coefficients (LECs) appearing above, such as p,p,K 
are related to the mass density, the shear modulus, and the compressibility of the solid respectively. 
They determine the velocities of the phonons in the solid phase. Similarly, the velocity of the phonon 
in the pure superfluid case is given by v^. In the presence of strong coupling between the solid and 
superfluid these coefficients are modified. For example, the coefficient p in Eq. [T] differs from the usual 
mass density of the pure lattice component due interactions that entrain the superfluid, and the mixing 
coefficient g m i X couples superfluid and lattice dynamics. As we will show Galilean invariancc relates 
Smix to the modifications of p and due to entrainment [9]. An analysis of these modifications in 
the context of the neutron star crust due to the underlying interaction between neutrons and protons 
was the original motivation for this study. In this case, the mixing coefficient g m i x is relevant for heat 
transport properties in the inner crust [TU], and the eigenmodes of the coupled superfluid-solid system 
could play a role in explaining the observed quasi-periodic oscillations in magnetars flares 

We will present a general proof that the functional form of the lowest-order Lagrangian is completely 
specified by the thermodynamic pressure in the presence of constant external fields that couple to the 

1 The U(l) symmetry is related to particle number conservation and we will refer to this as a phase symmetry. Its 
breaking simply refers to the choice of a ground state: total number is conserved and the continuity equation remains 
valid. 
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conserved densities and currents in the system. The derivatives of the pressure with respect to these 
external fields determine the low energy constants. Non-perturbative techniques such as Quantum 
Monte Carlo or Hartree-Fock techniques may be suited to calculate these thermodynamic functions. 
For example, the energy as a function of the density for a non-rclativistic uniform Fermi gas at unitarity 
was calculated using Quantum Monte Carlo techniques in [12] . (For a recent example of the calculation 
of the LECs in relativistic supcrfluids, see [IB]-) Since these derivatives of the thermodynamic functions 
are related to the long-wavelength limit of static correlation functions of currents and densities, their 
direct calculation using non-perturbative methods also provide the needed LECs. 

The outline of the paper is as follows. In Section [2] we outline the formalism and define notation. 
In Section [3] we revisit the effective theory of a neutron superfluid in the absence of any lattice, derived 
earlier in Rcf. |14[ I15j . This serves as a pedagogic warm-up before describing the more complicated 
system including the lattice. Here, we prove that lowest order lagrangian is determined by the ther- 
modynamic pressure. In Section 2] we derive results for the relevant case of the combined neutron and 
proton sectors. In Section [5] we focus on the applications of our formalism to the neutron star crust 
and use simple estimates for the mixing coefficient in the neutron star inner crust to determine the 
resulting eigenmodes of the longitudinal lattice and the superfluid phonons. Here we also comment 
on the connection with previous work on the elastic properties of LOFF phases |16j . We present our 
conclusions in Section [6] 

2 Symmetries of the underlying Hamiltonian 

The prototypical system we consider here is composed of two conserved species of particles. In the 
neutron stars, this would be the strongly coupled many-body system of non-relativistic neutrons and 
protons. In what follows we will continue to refer to these two components as neutron and protons but 
it should be understood that our considerations apply more generally. We represent the action for the 
system abstractly as S^n, tpp] where ^ n and *B P , are the neutron and proton fields. Even though the 
particles may be non-relativistic, we will work in a Lorentz covariant form and take the non-relativistic 
limit at the end of the calculation. We will denote the spatial indices by Latin characters (a, 6, c, i, j, 
k) running from 1 to 3 and the full space-time indices by Greek characters (fi, is, a, A) running from 
to 3. 

The theory describing the neutrons and protons is invariant with respect to global phase rotations 
of the neutron field, ^ n (x) — > exp(—i9 n )^! n (x). The conservation law associated with the symmetry is 
the conservation of neutron number. Similarly, independent phase rotations of the protons gives rise 
to the conservation of proton number. We note that protons and neutrons are separately conserved on 
timescales small compared to the weak interaction time r wea k — 1/(^f tib T 2 ) where tig and T are 
the baryon density and temperature of the system. For typical conditions the timescale for low energy 
dynamics is teft — 1/T and correspondingly the ratio t^ft /"Tweak — G 2 f ubT <C 1. We shall refer to 
the corresponding conserved currents as j£ and jp 1 respectively. The action S is also invariant under 
space-time translations and spatial rotations and the conserved current associated with translations is 
the stress energy tensor . 

To analyze the constraints provided by these symmetries on the low energy effective action, it is 
useful to introduce external fields that couple to the conserved currents (see e.g. discussion in Rcf. [17]). 
For the internal symmetries we add to the action S[*ff n , ^> p ] source terms of the form J d 4 x j£(x)A™(x) 
and / d 4 x j^(x)AP l (x). If the external fields arc allowed to transform appropriately under local U(l) 
transformations, this procedure promotes the global symmetries to local symmetries. This is equivalent 
to converting all partial derivatives <9 M to covariant derivatives — iA^. 

A similar extension of space-time symmetries is slightly more subtle. It requires extending space- 
time from a flat space-time to a curved space-time, and writing the action in a form that is general 
coordinate invariant. Indeed, the non-relativistic theory of neutrons and protons may be seen as the 
non-relativistic, flat-space limit, of a fully relativistic, general-coordinate-invariant action. The external 
field that couples to the stress-energy tensor is a deformation of the metric, 8g py . 

Therefore, we start with an action of the form, S[$/ n , ^> p , A^,A^, g^] which is invariant under general 
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coordinate transformations, 



^ X )^g^)=gP'( X ) — — (2) 



A p {x)^ A p {x')= A a {x) 



dx p dx" 
dx'v 



dx a 



these transformations include both rotations and boosts as special cases of local space and time trans- 
lations. The action is also invariant under local phase rotations of neutrons, 

* n (x) -> %(x) = exp(-i0„(x))*„(x) 

A-(x)^A';(x)=A-(x)-d,0 n (x), 

and local phase rotations of protons, 

%(x) -> %(x) = expH^CaO)*^) 

A*(x)^A'*(x)=A*(x)-d ll P>(x). 

The reason for extending the global symmetries to local symmetries is that correlation functions of 
the conserved currents can be analyzed very simply by taking appropriate functional derivatives of the 
generating functional W[A^, A^, g^} with respect to the external fields A™(x), A^(x) and g^(x). The 
generating functional is defined in the standard path integral representation as 

= Z[A^,AP,g^}; 

and Z[A™, A?, g^ v ] is thermodynamic partition function. For example, the derivative with respect to 
the zeroth component of the external field A^ defines the number density as given by 

5W[A™AZ,guv} 1 SZlA^APg.J 

In order to evaluate correlation functions in an equilibrium state with specified number density, the 
functional derivatives with respect to A^ (x) are evaluated at specific values corresponding to appropri- 
ate chemical potentials as required by Eq. [51 For neutrons A™(x) = A™ = (/x n + m n , 0), where fi„ is the 
usual non-relativistic chemical potential and similarly for protons A p {x) = A? = (fj, p + m p ,0) and \i v 
is the corresponding non-relativistic chemical potential. Moreover, functional derivatives with respect 
to g^vix) are evaluated at space-time metric g^vix) = g^ w . For the pure neutron sector (Section [3J|, it 
suffices to set the equilibrium metric to be the Minkowski metric, g^ u = r)^. In the case of coupled 
system, since the spatial components of a space-time independent metric specifies the lattice structure, 
we will allow g^ v to be more general in Section [5] In the following wc discuss specific cases in which the 
ground state |fi) spontaneously breaks number and translation symmetries of the underlying Hamilto- 
nian. First, in Section [3] we discuss the simple case of a supcrfluid which breaks the U(l) symmetry 
associated with number conservation, and subsequently in Section U we discuss the system of interest 
where both the global U(l) and space-time translation symmetries are simultaneously broken. 

3 One component superfluid 

3.1 Fields and the effective lagrangian 

To illustrate the main ideas we first consider a single component supcrfluid such as degenerate neutron 
matter where attractive interactions lead to the formation of Cooper pairs and a transition to a superfluid 
state. Here, the two-neutron operator has a non-zero expectation value and in equilibrium 

<Q|*„(x)vl£(x)|fi) = C 7 5 6(x) = C 7 5 |8| . (7) 
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Since the phase of the condensate changes on making a global phase rotation on the condensate 
spontaneously breaks U(l) n to Z^- The corresponding Goldstonc boson field <fi(x) is given by the phase 
fluctuations of the order parameter. 

The field 4>(x) transforms nonlincarly under U{l) n transformations with parameter 9 n : 

4>^^ + e n . (8) 

The effective Lagrangian for 4> can be in principle obtained by integrating out the heavy modes 
corresponding to the gapped fermionic excitations from the system in a Wilsonian approach. 

The partition function in presence of external gauge fields admits the following low-energy represen- 
tation: 

Z n [A™, 9fiV ] = f[d* n ]e iS ^> A >^ — > J [dcfie 15 *^^^ . (9) 

The symmetries of the underlying fundamental theory impose stringent constraints on the form of the 
effective lagrangian defined by S c s = J d 4 x s f^g £ c ff . Global U(l) n symmetry implies that the Goldstone 
boson can occur only through the derivative d^cj) and local U{l) n symmetry (Eq. [3]) implies that j4™ 
and 4> can appear in C e s only in the combination 

D^(x)=d^(x)+A^(x). (10) 

Since we are working in a covariant theory, C e g should transform as a scalar density under general 
coordinate transformations. Therefore, the effective lagrangian can only be constructed from building 
blocks like D^<fi, D V D^4> etc., with all indices contracted. 

In this work we use the power counting scheme proposed by Son and Wingate |15| to organize terms 
in £ e ff- We define the power of an operator as the difference between the number of d's and the power 
of <p. I.e. all terms of order (d) m (4>) n {m > n) with the same value of m — n are considered to have the 
same order, m — n. A™ has the same order as d^cj), i.e. order 0, and g^ v also has order 0. Therefore, 

C cff [D^g^]^C Q [D^D^}+C 2 [(D v D^) 2 ,. ..] + ... (11) 

and the leading order lagrangian £ is an arbitrary function of the building block X = g^ u D^D^cf)- In 
what follows we will focus on the lowest order lagrangian and relate its coefficients to the thermodynamic 
derivatives. 



3.2 Thermodynamic matching 

We now relate the functional dependence of C$(X) on X to the functional dependence of the thermo- 
dynamic pressure P(n n ) on the chemical potential fi n . We refer to this result as to "thermodynamic 
matching" . Although this result is not new [T4j [15] , here we provide a derivation that can be generalized 
to other, more complex patterns of symmetry breaking. In fact, we will use the generalization of this 
result in Section 0] when we consider the simultaneous breaking of translational and particle-number 
symmetry. 

We recall that the thermodynamic interpretation of the functionals Z[A^,g^] and WlA™, g^ v ] at 
constant external fields (A™(x) = A™ and g^ v {x) = g^ u = i^u) is given by the relation 

Z n [A^,ri^} = e tW - [A ^' n ^ ] = e - iVTQ " = J [d<}>} e iSm[D^4-,v^] y ( 12 ) 

where f2„ = (il\H — A™j^|f2) is the free energy density, V is the volume and T is the extent in the time 
direction. Here is the ground state of the Hamiltonian modified by the presence of the external source 
A™ (H is the Hamiltonian density). For the specific choice A™ = (m n + /i„, 0) and at zero temperature, 
\Q) is the many-body ground state at chemical potential [i n and fl n = (^„+m„)j" = — P(/x n ), 

where P{fJ. n ) is the usual thermodynamic pressure of the system. 

The low energy effective action S c s [D^(j>, r}^] is a function of the external fields and the Goldstonc 
fields. It contains all the quantum dynamics of the high energy Fermionic modes encoded as low energy 
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coefficients. In order to evaluate the partition function Z[A™, g^ v ] we expand the effective action about 
its saddle point 4>o which satisfies 

SScsiD^ix),!]^} dC e s(D tl <j>(x),ri l j U ,). . , 

5^5 ^ ~ — dd^) — Uo ' (13) 

minimizes the Euclidean action, and is well behaved at infinity. For general external fields A" (x), the 
solution to Eq. [T3]is a functional of the external field, ^>o[^]- However, for our homogeneous and static 
system with constant external fields the well behaved solution is 4>o — 0. Expanding about this point 
we can write 



1 



SeslD^r)^] = S c ffL=o + ^ J d 4 xd 4 x'ip{x)ip(x') Uo + ■■• , ( 14 ) 

where tp = (f> — 0Oi an d thus Eq. [12] can be evaluated as a loop expansion 

giWi-ioop = J [rfy] e '(5 / J'^VWy^') ^;.)^) Uo+-) ; (16) 

where we have explicitly displayed only the quadratic (Gaussian) part of the functional integral in [dtp] 
which corresponds to the one-loop approximation. Let us now discuss this loop expansion in light of the 
EFT power counting. The key observation, which is a generic feature of low-energy effective theories [18] , 
is the following: within the momentum (gradient) expansion of the EFT, loop diagrams generated by 
Co are higher order than tree- level diagrams with vertices from Co- in our case, one-loop contributions 
to phonon amplitudes are suppressed by four powers of momenta compared to the tree graphs generated 
by £q as shown by Son and Wingate |15] . Using the above considerations we can write 



W[A»] = I d 4 xC eff ((D^o[^})^^) + Wx_ Ioop (^) + ... 

(17) 

Co(Xo) + C 2 [A%) + C± [A™] + Wx-ioop(^) + ... , 



d 4 x 

where X = D^oD^^o- Cq{Xo) in Eq. [T7]is the leading term (O(p )), the second term is of 0(p 2 ), 
the third and fourth are 0(p A ). The contribution of O(p ) involves either no derivatives on the external 
fields or two derivatives compensated by a Goldstone propagator of 0(p~ 2 ). Higher order contributions 
necessarily involve derivatives acting on the external field A"(x). So we arrive at a very important result: 
for very long wavelength external field (A™{x) — > constant), only the first term in Eq. [T7] survives., i.e., 

W[A"] = [ d 4 xCo{AlA^ n ) = VTCo{A n A^ n ) (18) 



Now recall that A^A^ = (m„ + p n ) 2 so that fi n = y/ A^A^ — m n . Moreover, Co depends on D^cj) only 
through X. At the classical solution, for constant external field, one has X Xo = A^A^. So we have 
from Eqs.nUandCEl 



Co{Xo) = PtyAnA* - m„) = P(VX - m„ = fi n ) . (19) 

The above relation fixes the functional dependence of Co on the variable X once the functional depen- 
dence of P on p n is known. So in general we have: 

Co(X) =P(Y = VX-m n ) . (20) 



Finally we note that in the non-relativistic limit the relevant building block takes the form Y 

2m, 



V(m n + fi n + d <f>) 2 - (dit) 2 - m n ~n n + do<t>- HUGHES]. 
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3.3 Identifying the low-energy constants 



Eq. [20J gives us the complete expression of the superfluid lagrangian to the lowest order in (m — n) . 
Expanding the function £o(4>) in powers of the Goldstone fields, one can read off the phonon kinetic 
term and self-interaction vertices: 



£eff[4>] = A) [(/in 
1 d 3 P 



m n y 



dP 

dY 



3! d ya (^) 3 + » 



(do® 



1 



1 d 2 P 
V.dY 2 
d 2 P 



(do® 2 - 



1 



dP 



2\(p n + m n ) dY 
1 dP 



Qi„ + m n ) dY 2 (ji n + m n ) 2 dY 



(21) 



(do^d^di® + 



where all derivatives are evaluated at Y = fj,„ . The above expansion makes it clear that the low energy 
constants of the theory are then given by the derivatives of the pressure with respect to the chemical 
potential. Eq. [5T] also serves as a starting point for a fluid dynamical study of the system once one 
realizes that —di4>/(m n + fj, n ) is the velocity of the superfluid. (See Ref. [2D] and references therein.) 

Alternatively, one can show that the thermodynamic derivatives are related to the static correlation 
functions involving the neutron charge and current operators. For example, 



dWn 

9A« 



dP 



(22) 



Similarly, 



d 2 W 

dA n dA n ]ieq ^ i I d4x ( T { J n( x ) J n(°)})oonnected 



dP, 



(fin +m n ) dY 



H n + m n dY + dY 2 ) leqV 11 



(23) 



A more careful analysis taking A to be long wavelength but not quite constant shows that the current- 
current correlation function thus obtained is the momentum independent part of the transverse-current — 
transverse-current correlation function. 

Finally, in the non-relativistic limit, separating the space and time components we obtain the well 
known results 



d 2 W n 

d 2 W n 
dAfdA? leq 

1 J 



d 2 P 
dY 2 
1 dP 
m„ dY 



(24) 



4 The superfluid and lattice phonon lagrangian 

4.1 Fields and the effective lagrangian 

A crystalline ground state of the system spontaneously breaks translations and rotations. The proton 
number density acts as the relevant order parameter. Since space-time dependent translations include 
rotations as special cases, it will suffice to consider the breaking of the abelian group G of spatial 
translations to the subgroup H = {TA containing discrete translations by multiples of lattice basis 
vectors, b. The generators of G are the components of the three-momentum P a given by space integrals 
of the energy-momentum tensor components T 0a (x), and momentum conservation takes the local form 

d^ a (x) = 0. (25) 

The Goldstone effective fields can be chosen as space-time dependent coordinates £ a ( x ) ( a = 1) 2, 3) 
of the coset space G/H: 

1 (x)=e* a{x)pa , 7 e G/H . (26) 
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The nonlinear action of the translations group on the Goldstonc fields is specified by: 



Xb -> x'b = x b + a b , £b(x) -> ('b(x') = 6>(z) + a b , (27) 

as can be verified by left multiplication of 7 £ G/H with g = e l ° hph G G. 

Promoting the global symmetry to local symmetry, a generally covariant formulation [21[ 114] of the 
phonon dynamics in the background metric g^ u can be readily achieved by introducing a set of fields 
that transform as scalar fields under the general coordinate transformations of Eqs. [2] 

z a (x) =x a -C(x) a = 1,2, 3. (28) 

The fields z a (x^) can be thought of as one particular choice of body- fixed coordinates d of a material 
point located at a; M = (t,x). With this choice the body-fixed coordinates coincide with the "laboratory" 
coordinates (x a ,g a b) when the displacement field £ a vanishes. 

As in the pure neutron case, the partition function in the presence of generic external fields A™(x), 
A v Ax), g^(x), admits a low-energy representation in terms of the four Goldstonc modes <j) and £ a : 

Z[A^,Al,g^]= [[d^njld^p]^ 3 ^' 91 '^'^' 9 ^ -> [[d(t>}[d^ a ]e iSs!t[ ' t> ' e ' A ^ A >''" ] . (29) 



At the end, we will evaluate the partition function for space-time independent external fields .A™, A^ 
and g^, specifying a particular density, and lattice shape for the system. 

iS e ff represents the effective action of £ a (or equivalently z a ) and </> in the presence of external fields. 
We can organize the terms in S c s according to the same power counting introduced earlier in our 
discussion of the supcrfluid, i.e. in increasing difference between the number of derivative operators and 
the Goldstone fields, 



SafaCA^A^g^] = J d^x^g-yCoid^d^z^A^A^g^) (30) 

+ C 1 {D v d^D v d^z a ,D tl Al. ..) + ...] . (31) 

Symmetries impose powerful constraints on the form of £ c ff- Since z a transform as scalars, d^z a 
transforms as a contravariant vector. The building blocks of the scalar function C e s are scalar com- 
binations of A"(x), A p {x), g^ v (x), dfj,(j), d^z" 1 , and their covariant derivatives. Symmetry under phase 
rotations of the neutrons, Eq. [3J implies that A™(x) should appear in a combination such that the trans- 
formation A™{x) — > A™(x) + d^9 n {x) leaves the effective action invariant. The same is required for the 
protons. In the pure neutron case we found that gauge symmetries implied that d^cj) and A™(x) could 
appear only in the combination D^(f>. To lowest order in the power counting, the scalar combinations 
that can be constructed from the gauge invariant combinations are 

X = gTD^Dvtj, (32) 
W a = g^D^d v z a (33) 
H ab = gT 8 fl z a d v z b . (34) 

In addition to these building blocks, other possibilities arise in the mixed case that were not present in 
the case of a pure neutron supcrfluid. The following terms only change by a total derivative on making 



2 These are the coordinates in a frame frozen in the body of the solid. If one follows a material point in the solid, its 
coordinates in this frame remain constant. 
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gauge transformations defined in Eqs. [31 0] 

J d^e^e^fc, A*(x) + C 2 d^{xj) {d u z a {x)d a z b {x)d x z c {x)) . (38) 

Hence, the most general form of Co is 

C (d^,d„z a ,A^Al, 9llv ) = f(X,W a ,H ab ) 

+ -^=e^ aX e abc (d Al + C 2 d^) (d,z a d <J z b d x z c ) . ^ 

The term proportional to C 2 is a total derivative that becomes relevant only in presence of non-trivial 
topological configurations (vortices) for the field <f> [5] . This would be relevant for calculations of vortex- 
phonon interactions but for now on we disregard this term and restrict our discussion to vortex free 
configurations. 



4.2 Thermodynamic matching 

Extending the analogy with the neutron superfluid case further we can relate Co to the free energy of 
the neutron-proton system. The free energy ^l[A^, A^, g^] is proportional to the log of the partition 
function. Following the discussion in the pure neutron case, one can show that for space-time indepen- 
dent external fields, g^ix) = g^, A^{x) = A?, A" (x) = A™, f2 it is also equal to Co evaluated at the 
classical solution <fi\o = 0, £ a |o = 0. Hence, 

ZlA^A*,^] = = e-ivm^M = ei m D (oi;i;J S *„) > (40) 

where VT — J d 4 Xy/—g. For this choice of the many body ground state, Xo = A p \A^ ,n 1 Wq = and 
H ab = gab Therefore, 

-n\pL n ,v P , g^} = f(X = Xo, W a = 0, H ab = g ab ) + Ci (jh + m p ) . (41) 

The constant C\ can be determined from the requirement that J^- = -§== . Thus, we see that C\ 

is the density of protons for a configuration whose metric has determinant —1. Symbolically, C\ = rip 1 , 
where rj^ is a particularly convenient choice (also see footnote 5) for a metric with determinant — 1. 

When we consider the functional form of /, we encounter a feature different from the previous 
case where we considered the pure neutron superfluid. There, we were able to determine the complete 
dependence of the function / on its arguments from the free energy function fi (Co(X) = P(Y) = 
— Q n (Y)), i.e. from a calculation of the partition function with the specific form j4" = (m n + fj, n , 0) for 
the external field. In the mixed case, however, since D M </><9' i z a \ eq = for ^4™ = (m„ + /i n , 0), it is not 
possible to calculate the dependence of / on W a from the free energy calculation in this external field. 
To determine the dependence of / on W a one needs to evaluate the partition function Z for a space-time 



The term 

d A x < ^e^< lX 6 abc (A™{x)YdvZ a (x)d <T z b {x)d x z c {x)) (35) 



3! 

can be rewritten as, 



J d 4 x^ [3!^ 9 det(H) - e^" x e abc (d ll 4>(x))(d v z a (x)d l7 z b (x)d x z c (x))\ , (36) 



where, 

X W aT 



JH 



W a H 



a b 



(37) 



This shows that any term proportional to C3 can be reabsorbed by a redefinition of the function / and the coefficient C2 ■ 
4 This fact is intuitively understandable. In the non-relativistic limit 1141 we have W a ~ rn n (— ^— d a <t> — &o£ a + 
■^—9i4>di^ a ) = m„(v® — do^ a — ^ n -V£ a ) which is the relative velocity between the neutron superfluid and the proton 
clusters. The dependence on W a therefore represents the interaction between the superfluid neutrons and the lattice when 
they are moving relative to each other, and can not be calculated by a ground state evaluation of the free energy. 
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independent external gauge field A™(x) that has non-zero spatial components, A™ = {/J, n +m n , Aj). This 
gives D^o = A™, X a = A™M l n and W§ = g ah 'i™ = A a ". Then, 

-iilA^A*,^] = £0(0,^,^,4,5^) (42) 

= f(X = X ,W a = W a ,H ab =g ab ) + -L=n;(fi p + m p ) (43) 

By calculating the free energy for various A™ and g ab we can map out the functional dependence of / 
on X, W a and H ab . Finally, noting that Q = (Q\H - A^j£ - n p {fx p + m p )\Q) (H is the hamiltonian 
density), the term proportional to n p cancels out from both sides and the function / is given by, 

f(X , W a , g ab ) = (fi|J&tf - H\Q) . (44) 
The generalization of Eq. 00] with the full A n is simply 

Z[A n ,A p ,g^] = e iWlA ?'- A *' S ^ ] = e-iVTUlA^^g^] = jVTC o (0,SZ,A»,A*,g hv ) _ (45) 



4.3 Identifying the low-energy constants 

Expanding the function Co in powers of the Goldstonc fields <p an d £ a , one can read off the phonon 
kinetic term (including kinetic mixing among the £ and <f>) and self-interaction vertices. The expansion 
in the 4> field can be done as in Section [3] The expansion in £ a is performed about the undeformed 
equilibrium configuration with g^ v = t]^ and £ a = 0. Deviations from the equilibrium shape are then 



signified by £ a 7^ 0, which gives H ab = 
equilibrium, X = X and W a = 



n 



ah 



AH ab with AH 



ah 



-(d a £ b + d b £ a ) + dn£ a d v Z b if. At 



0, and deviations from equilibrium are given by X — Xq 



AX = 2(/x„ + m n )d <j) + d^d^if" and W a 



W§ 



AW a 



To second order in the fields, the expansion of / is 
f(X, W a ,H ab ) = f(X ,0,v ab ) 

dX 
1 



+ 



(2(/i, 

eq 

d 2 f 



2 dW a dW b 
df 



QJJab 
d 2 f 



eq 

AH ab 



(-( Mn + m n )d°C + 



1 a 2 / 



1 



2 dH ab dH cd 



2! dX 2 
AH ab AH cd 



(2(/i„ + m„)aV) 5 



+ m n )d°e + d 



(46) 



dH ab dX 



AH ab {2{^ n + m„)9V) 



We have simplified the expansion above by taking J^ a le( 



0, 



and a„.iL 



0. 



This would be the case for any crystal with reflection symmetry, for example a cubic crystal. For a 
cubic crystal, one can further simplify the expressions by using symmetry under rotation by ^ along 

the axes This trives ^ f — 1 32 / xab df _ 1 df tab j d 2 f _ 1 d 2 f xab W W P 

me axes, ims gives dwagw t - 33rar o , - 39H ^o , ana gxdHab - 3d H<"=dx° ' wncrc 
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-qjjtx + Qffuu - Finally, we make the non-relativistic approximation fj, n -C m n , fi p <C m p , 



and separate the space and time components 
With all these simplifications, 



f(X,W a ,H ab ) = f(X ,0,r ] 



ah \ 



2 df 



d 2 f 



3dH 
1 df 



ZdH c 



n idW c dW c _ 

1 d 2 f 



m n dY 

d 2 f 



d 2 f 



2 8H ah dH 



^(d a e+d b c)m d +d d n + 



3dW c dW c 
d 2 f 



n 3dW c dW c 



(47) 
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One key consequence of Eq. 27J is that the low-energy constants are related to derivatives of the 
function / with respect to X, W a , and H ah evaluated at the "equilibrium" point X = Xq, W a = 
Wq, H ab = 77 . In turn, due to the thermodynamic matching relation Eq.03] the low-energy constants 
can be expressed in terms of derivatives of the generating functional W[A™, A p , g^]. The analysis 
proceeds along parallel lines to the pure neutron case, but it contains a number of novel features, which 
we discuss in some detail below. 



4.3.1 Thermodynamic derivatives 

The first order derivatives of the functional W[A™, AP^g^ specify the number density of particles in, 
and the stress energy tensor of, the system: 

1 dW . . . df 



VTdA^ eq ~ v n/ ~ dY 



-2 dW 00 , J/ 



(T 00 ) = (m n + Ai „)^-/ (48) 



■ 2dw - -w- 2 y4-/K 



VTdg ab q 1 1 3dH c 

In particular, 2 jr 8 d J ac = — \(T£) — f. The second order derivatives of W have the following form, 



' 3 dH"<= 3 

1 d 2 W d 2 f 



- F A Ao 



VTdA%dA% ]eq dY 2 

1 d 2 W f J_ 9/ 1 d 2 f 1 „ = _ 

VTdA n a dA^ eq l m n dY 3dW c dW ci ' AaA " 



d 2 W d 2 f 



2 — FA g 00 



VTdA%dg 00 ,q dY 2 
-2 d 2 W , _ r df o l d 2 f 



1-—L - 9- J ]n ab = — TP a 



VTdA%dg ab ,q 1 dY 3dH cc dY i 
-2 d 2 W f 9/_ d 2 f 1 ab _ 

VT dA^dg bo leq ~ [ 8Y n 3dW *dW c 1 V ~ Aa9b0 (49) 

4 d 2 W , df_ 2 fl 
VT dgoodgoo Uq _ " d^ 2 ~ 9 ° 0ff0 ° 



VT dg 00 dg ab leq ~ [ 1 n dY + 3dH°c 3 "dYdH~ lV ~ 9oo9ab 

4 d 2 W , , £?/_oI_^l_ 2 g 2 / 1 q b _ j? 

VT dg a0 dg bo leq ~ 1 I+mn 8Y 3 dH cc mn 3dW c dW clV ~ 9a09bQ 



VTdg ab dg cd ,q w 3dH ee ' yi ' ' 1 111 dH ab dH cd 

= — TP 

— 1 QabQcd 1 

where all derivatives of / are evaluated at equilibrium. The second order correlations are proportional 
to the momentum independent part of appropriate time ordered correlation functions of the neutron 
current and the total stress energy tensor. Since they can be found simply from Eq. [49] by noting that 
j M and T^ v are obtained as the partial derivatives of W with respect to A^ and g^, respectively, we 
don't explicitly include them here. The expressions in Eqs. l46l and l49l look complicated but have simple 
physical interpretations, as we discuss below. 

4.3.2 The entrainment coefficient 

Here we note that the current-current correlation function for neutrons, g^„^„ \ eq is not simply pro- 
portional to the total neutron density but is instead proportional to — q^cq W c ■ We conjecture 
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that n b = g W ag W c > and this represents the number density of neutrons "bound" or "entrained" 

on the nuclei, ^p- — m n | q\y<^ow c 1S ^ nen interpreted as the number of "unbound" neutrons in the 
system. Indeed, from the coefficient of {di<j)) 2 in Eq. 27] we see that the current of the supcrfluid mode is 
proportional to nj — (n n — n ) , where nf is neutron density that can participate in supcrfluid transport. 

It is also reassuring that in Eq.|47]thc effective mass density of the "proton" clusters involved in lattice 
vibrations is correspondingly increased by m n rn,. This is easily seen by noting that the coefficient of the 

kinetic term is given by ^[f g^/cc + m n^ ow c dW a ^ = itf gjjee + fn n n^[. The picture of the inner crust as 
periodic clusters of ions and neutrons "entrained" on the clusters has been discussed previously [2] [3] . 
Our formalism confirms this intuition, and provides a field theoretic derivation of the entrained neutron 
density in terms of generalized thermodynamic derivatives. 



4.3.3 Relating the LECs to the stress and elastic tensors 

From Eq. 06] and the last relation in Eq. [49] one sees that the part of the effective lagrangian quadratic 
in gradients of £ a can be expressed in terms of gfj^g d leg- This result establishes a non-trivial relation 
between the LECs appearing in the phonon quadratic lagrangian and the stress tensor correlator that 
can be calculated in the underlying theory using non-perturbative methods. We can go one step further 
and relate the LECs to first and second order thermodynamic derivatives of the free energy with respect 
to the strain tensor (i.e. the stress and elastic tensors). This step relies on the relationship between 
the external metric <7 M „ and the strain of the crystal structure, i.e. its shape. We will find that since 
the strain has pieces both linear and quadratic in the displacement fields (see Eq. \ST\ below) the elastic 
constants are linear combinations of first and second order derivatives of the free energy with respect 
to the strain. 

Let us first recall a few basic definitions from the theory of elasticity. The elastic constants can be 
defined through thermodynamic derivatives of the free-energy (or internal energy) per unit mass with 
respect to the strain tensor s ab associated with deformations around some reference point: 

F = Fq — t ab s ab + X - X abcd s abScd + ... (50) 

In the above relation t ab is the stress tensor associated with the reference configuration, which we will 
take to be the equilibrium configuration. x abcd is known as the elastic tensor. For an equilibrium 
configuration in absence of external forces, t ab = 0, and the components of X abcd are simply the elastic 
constants. For an equilibrium configuration in the presence of external forces (t ab ^ 0), for example a 
solid under pressure, the elastic constants are linear combinations of t ab and X abcd . 
The strain tensor is defined in terms of the displacement fields ( a (x) as follows: 



1 fdC, a 


d( b 


dC dC 


2 \dx b 


dx a 


dx a dx b 



Sab 2\dx b + dx- + dx a dx b ) ' (51) 

The strain tensor has a simple geometric interpretation. Imagine choosing the body-fixed coordinates x a 
so that they coincide with the Euclidean (flat) laboratory coordinates when the body is in equilibrium. 
After a deformation specified by the displacement fields ( a (x), the body- fixed coordinate system will 
have a non-trivial three-dimensional metric, whose deviation from flat metric is specified by s a b- 

Vab = S a b -> 9ab = Vab - 2 S ab . (52) 

We now state the results that ensure the connection with the elastic constants, relegating their proof 
to Appendix[3] The main point is that the energy density f2[A n , A p , g] calculated using the path integral 
(Eqs. I29land l40l) in the presence of a space-time independent metric of block form 



9^ 



1 

9ab 



(53) 



5 Equivalently, one defines the thermodynamic quantities per unit volume of the undeformed 1221 system. The free 
energy per unit flat space volume element is given by F = \/—gQ = — W[A™, A^,g^^]. 
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is equal to the flat-space energy density in the lowest energy state |f2 s ) subject to the "deformation 
condition" (n g \£ a (x)\fl g ) = with Q(x) related to g ab by 



Jab = Vab - 2s a b(Cff) 



(54) 



and s a b(C) given in Eq. 1511 This result establishes a correspondence between the ground state in presence 
of (job 7^ Vab and a deformed configuration around the "true ground state" in the absence of external 
gravitational field, g ab = r\ ab . Therefore, by varying the external metric g ab we probe different deformed 
configurations of the system, with strain tensor related to g ab by Eq. 1541 



Identifying the free energy F[s] per unit volume in fiat space with F 
we have: 



VT 



bed 



X a 

For a cubic crystal 

-^abed 



S 2 F 



5Sab5s cd 



8 2 yf=g~n 



SSabSs 



cd 



S 2 W 



VT 5g ab 5g, 



cd 



■W[A™APg^] 



(55) 



_ Zji ) s ab S cd + p. (S ac S bd + S ad S bc ) + aS aba 



(56) 



where the term proportional to a is non-zero only if a = b = c = d and represents the anisotropy in the 
clastic coefficients. Now we have all the pieces to write the LECs in terms of thermodynamic derivatives. 
For convenience, we define P = — |(T^), and E = (T 00 ). 

4.3.4 The quadratic phonon lagrangian 

Making the identifications discussed above, we can write the quadratic lagrangian in a rather compact 
form. Neglecting constant terms and total derivatives, and using integration by parts to simplify some 
terms, we find: 



A. = l[-F AoAo ](d 4>) 2 - l - 



i 



-^{ F A a A b )Vab 



P + E+^L(F AaAb ) Vab 



m n F AaAb )r] ab 



(57) 



This form allows us to express the low-energy constants appearing in Eq. [T]in terms of thermodynamic 
functions and derivatives as follows, 



P = P + E+^(F AaAb )i lab 

K = K+-P 

H = p,- P 



f(j> 

*?mix — 



- F A A 

--^F AaAb i] ab 



(58) 



-^J-^{ F A a g ab + m n F AaAb )rj ab \ , 
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where 



K = (-^S abcd + U ab S cd + l5 ac 6 bd )X abcd (59) 
18 6 9 

P = {-ls abcd + h ac S bd )X abcd (60) 
6 6 

a = (ls abcd - U ab 5 cd - U ac 5 bd )X abcd , (61) 
6 6 6 

with X abcd given in Eqs. E3 EH 

The pressure term in the definition of the bulk (K) and shear modulus (/i) in Eq. [55] should not be 
cause for concern. Its origin can be traced back to the term linear in s ab in Eq. [50] which is present in 
a system at finite pressure [H] . 

This dependence of the elastic constants on the linear term in the Taylor expansion of the free energy 
with respect to the strain tensor s ab appears to be counter-intuitive. However, the key point here is 
that the strain tensor associated with a deformation £ a has parts both linear and quadratic in d b ^ a 
(Eq. [51]). Using this, one can show that the elastic constants are entirely determined by the quadratic 
terms in the expansion of the free energy with respect to the displacement field £ a (that appears in the 
combinations £ ab and d a £ a ). 



5 Applications 

5.1 Neutron star inner crust 

Here we apply the formalism to the inner crust of neutron stars and illustrate the importance of entrain- 
ment and kinetic mixing induced by the neutron-proton interactions. We revise the calculation of the 
mixing constant <? m i x (Eq. [1]) in Ref. |10| including the effects due to entrainment. In this earlier work, a 
somewhat arbitrary distinction was made between neutrons bound in the nuclei and the neutrons "out- 
side" . The interaction between the nuclei and unbound neutrons was modeled by a short-range potential 
V(r) = —2ira n i 6 3 (r)/m n where a n i was the effective neutron-nucleus scattering length. Here, using 
the results of the previous section we show that neither of these ad hoc assumptions are necessary as the 
LECs of the effective theory are simply related to generalized thermodynamic derivatives evaluated in 
the non-perturbative ground state. A first principles calculation of the LECs would require a numerical 
non-pcrturbative calculation. Such a calculation is beyond the scope of this study. In what follows we 
will use simple estimates based on earlier calculations to draw some qualitative conclusions about the 
role of interactions between the solid and the superfluid in the neutron star crust. 

From the preceding discussions the lowest order effective lagrangian for longitudinal modes with 
canonically normalized fields <f> = f^<j) an( i ti = y/p£i m the inner crust can be written as 

C = i(5 o 0) 2 - \v% (d^f + i(do6) 2 - \ vf(diii) 2 + <w d o d % l , (62) 
where the LECs defined in Eq. [58] can be written as 

j2 _ nf 2 _ K + (4/3)// _ 1 (F Aogab + m n F AaAb ) ru , , 

V 4> ~ 72 ' V l ~ ' anQ 9mix ~ o / = : ' ID-Jj 

m « f£ P 3 y/-F Ao A P 

Here we note that (—F AaAb r] ab /3)f/m n = (n n — n b )/m n where nt is the density of "free" neutrons that 
participate in superfluid motion, n b is the number density of neutrons entrained by the lattice and n n 
is the total neutron number density. Further, since neutrons and ions remain non-relativistic in the 
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neutron star crust, the LECs simplify 



p= E + P+^(F AaAb ) Vab 




Non-Rel , 



(64) 



Non-Rel, 




d 2 f 



1 




dn. 



dn p dp, 



n 



where the hybrid free energy function is 



f(^n,n p ) = {fJL n + m n )n n (p n ) 



E(fj, n ,n p ) . 



(65) 



Here n p is the proton density, K and p are respectively the bulk and shear moduli of the combined 
system, and we have ignored the small contribution due to the LEC a that encodes the anisotropic 
contribution. In |10j a simple estimate of g m i x was derived but failed to include the contribution due to 
entrainment effects. We have verified that the result in Ref. [10] can be recovered by setting rib = in 



In Ref. [10] , assuming that the effective interaction between the unbound neutrons and ions is weak, 
it was found that /? = m ri kpj'K^ and in = rif / (rn n f^), where fcp and rif = k F /3ir 2 are the Fermi 
momentum and number density of unbound neutrons. The speed of longitudinal lattice vibrations 
was approximated as the Bohm-Staver sound speed. The longitudinal sound speed is given by vi = 
\J Ki/p where Ki — p(d(Pi e )/dp) is the bulk-modulus of the electron-ion system. To calculate the 
longitudinal speed in the Bohm-Staver approximation, the total pressure of the electron-ion system 
is (well) approximated by the electron pressure P e , and the mass density of lattice is taken to be 
p = m n A where A is the number of bound nuclcons in the ion. Interactions between nucleons will modify 
these simple estimates quantitatively. Qualitatively, the effect of strong neutron-proton interactions is 
the induced mixing between longitudinal lattice phonons and the supcrfluid modes. This interaction 
is characterized by the dimensionless LEC, g m i x , which in turn depends on two contributions, one 
proportional to n p (dn n /dn p ) and the other proportional to the entrainment parameter rib- 

In the neutron star context, both of these quantities can be calculated using phenomenological 
models. The ground state structure which specifies the profiles of nucleons is obtained by solving 
the single-particle equations in the Wigner-Seitz (WS) approximation [I] or more realistic boundary 
conditions that reflect the cubic lattice structure J25[ [26]. For a given volume of the WS unit cell 
VvvS) these calculations determine the number of bound neutrons (Nb), protons (Z), and the total 
number of neutrons in the cell (-/Vws)- They also determine how p n and p p vary with neutron and 
proton densities. The first contribution to g m i K is found by noting that n p (dn n /dn p ) = n p f 2 V np where 
V np = (dp n /dn p ) is the effective interaction between neutrons and protons. The other contribution 
is related to the density of bound neutrons and a naive estimate would suggest that n b = Nb/Vyys- 
However, as discussed earlier in section 14.3.21 and in Ref. [9], the number density of neutrons that 
effectively move with the nucleus is defined through the static limit of the current-current correlation 
function n = —FA a A b Vab/3 = (n n — rib)/m. This correlation function has been computed earlier for 
neutrons in the background of a static periodic potential designed to mimic the neutron star crust in 
Refs. [2] [25]. In these calculations k = (n n — n*)/m* is defined in terms of an ad hoc but convenient 
quantity called the effective mass m* of unbound neutrons, and the average number density of neutrons 
with energy greater than zero is denoted by n* . Thus the LEC 



where n n = -/Vws/^ws is the average neutron density in the cell. 

We now turn to a simple illustration of how mixing affects the propagation of longitudinal modes 
in the crust. For this purpose it would be ideal to compute the three LECs («</>, fj, 3mix) from a self- 
consistent underlying microscopic model using Eq. 1631 However, such a calculation is beyond the scope 



Eq.[Ml 



rib = n n - 




(66) 
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Figure 1: The velocities of the two eigenmodes. The dotted lines are vi and ignoring mixing. 



of this work and we adopt a less rigorous approach where we use the results of Ref. [lOj for the velocities 
of the superfluid and lattice modes in the uncoupled system, and assume that rib ^> n p (dn n /dn p ) and 
m* /m ~ 1. Simple estimates support our expectation that the dominant contribution to <? m i x is due to 
rib- In this case, 

5mix ^ «0 7= = (67) 

\J{rib + n p )nf 

Our second assumption m* /m ~ 1 is likely to be invalid in some regions of the crust |27j . Nonetheless 
to simply illustrate the role of mixing we have set m* = m and plan to return to a fully self-consistent 
calculation in future work. 

In terms of the canonically normalized fields the kinetic terms in Fourier space has the form, 

S = I E0(-*) U{-k)) ( kl ~ft] !? xk tl ) ( ffn ) (68) 
2^ V 5mixfc |k| kl -«, a k 2 J V k.£(fc) J v ' 

The velocities of the two eigenmodes can be obtained by diagonalizing the matrix in Eq. HU The results 
are shown in Fig.Q]where the solid curves incorporate mixing effects due to a finite g m i x given by Eq. 1671 
and the dotted curves show the uncoupled case with g m - lyi = 0. 

In contrast, the speed of the transverse lattice modes are unaffected by mixing and is given by 

(69) 

Here, only entrainment effects play a role in the propagation of transverse lattice phonons, as was 
previously pointed out in Ref. [5]. 




5.2 Crystalline superfluids or LOFF-like phases 

Other systems of phcnomenological interest where this low energy theory applies are the LOFF phases [5j 
[6]. Here, attractive interactions between two species of fermions leads to pairing at the Fermi surface 
but with a pair condensate (1P11P2) which is spatially inhomogcncous in the ground state. A mismatch 
in the Fermi-momenta of the two interacting species in the absence of pairing, disfavors the formation 
of zero-momentum Cooper pairs and instead pairs with finite total momentum are favored. These pairs 
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condense to form a ground state that breaks translation symmetry and can be written as a sum over 
plane- waves, 

(Mx)M*)) ~ A(r) = A > ( 7 °) 

{q°} 

where A is the gap parameter. The magnitude of the momentum |q a | ~ Skp where SUf is the splitting 
between the Fermi momenta of the interacting species. The magnitudes of the momenta and their 
spatial orientation is determined by minimizing the total free energy and this set of momenta specifics 
a crystalline ground state. The LOFF phases can in principle be realized in cold atomic Fermi gases 
where a splitting between Fermi levels can be achieved through a population imbalance [28] and in dense 
quark matter where a Fermi level splitting arises naturally [29] 130] . 

In dense quark matter, pairing between different flavors of quarks can play a role in determining 
the ground state structure. The relatively large strange quark mass, and charge neutrality induce a 
splitting between the Fermi energies of up, down and strange quarks. The expected splitting between 
the Fermi energies is 5p ~ m^/(8fj, q ) where m s is the strange quark mass and [i q is the quark chemical 
potential. At moderate densities, such as those realized in the neutron star core where fi q ~ 400 MeV, 
this splitting between Fermi energies can favor LOFF phases in quark matter with spatially varying di- 
quark condensates with a crystalline structure when <5/i ~ Ao/y/2 [3TJ GUI [2H] > where Ao is the gap in the 
absence of Fermi surface splitting. As this ground state breaks the same symmetries as those discussed 
in Section [2] these phases are amenable to the same low energy energy effective theory formulation. 

Several aspects of the low energy theory of crystalline phases in quark matter have already been 
described in Ref. [34] . In Ref. [16], the coefficients of the "lattice only" ((f) = in Eq. [TJ effective theory 
were computed microscopically in a Ginzburg-Landau expansion [33] . This work was primarily focussed 
on the shear modes and showed that both the kinetic coefficient, p, and the elastic constants [i etc., 
were of the order p q A 2 where A is the pairing gap parameter. The mixing between the longitudinal 
lattice phonon mode and the supcrfluid mode was mentioned but the relevant mixing coefficient was 
not calculated. Using the same techniques, we have estimated the mixing coefficient and we find that 
in the regime where LOFF-like phases are favored 

loff _ A (71) 

ymix c v 1 / 

0/1 

Similarly, the coefficients for the "superfluid only" (£ a = 0) sector can also be computed (some were 
calculated in [35]). For A <C Hq we expect the coefficient /| ~ fi q , corresponding to the density of states 
near the Fermi surface for a rclativistic system. Our simple estimates here show that strong mixing 
between the supcrfluid and the longitudinal mode can be realized, with important implications for hy- 
drodynamic oscillations both in the context of dense quark matter and trapped imbalanced Fermi gases, 
where LOFF phases may also be potentially realized. Definitive results require a rigorous derivation of 
the low-energy constants and such a calculation is being pursued and will be reported elsewhere. 



6 Conclusions 

We have studied a low energy effective theory describing phases of matter that simultaneously break 
translational symmetry and a number conservation symmetry. U(l) phase invariance and general coor- 
dinate invariance restrict the combinations of terms that can appear in the effective lagrangian. We have 
shown that the lowest order lagrangian (featuring equal number of derivatives and Goldstone fields) is 
determined by the derivatives of the thermodynamic pressure with respect to the external fields such 
as the chemical potential. While this was known in the case of one supcrfluid system [14[ 115] , here we 
have provided a different proof for superfluids and we have generalized it to the mixed system. The two 
main results of this paper, Eqs. [49] and 1571 provide a useful framework for computing the low energy 
dynamics. 

Our thermodynamic matching relates the LECs to thermodynamic derivatives of the free energy 
with respect to external fields (chemical potential, vector potential, background metric). We have also 
pointed out the relation between LECs and correlators of the U(l) current and the stress tensor at small 
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momenta. Both approaches might be pursued in future non-pcrturbative calculations using many-body 
techniques such as Skyrme Hatree-Fock and Quantum Monte Carlo. 

As a concrete example of phenomenological interest, we have considered matter in the inner crust 
of a neutron star, updating a previous estimate of the parameter characterizing the kinetic mixing of 
supcrfluid and lattice phonons. We also discussed briefly how this formulation would apply to the 
crystalline superfluids or LOFF-likc phases and highlighted the role of mixing between the modes in 
these systems. 

Finally, we note that the formalism that we have set up here can be applied to study the low-energy 
dynamics of other physical systems with several spontaneously broken symmetries, such as a system 
composed of two superfluid species. 
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A Q[A n ,A p ,g] and energy density of deformed states 

In this appendix we show that the energy density il[A n ,A p ,g] = — W[A n , A p , g]/(VT) calculated using 
the path integral (Eqs. [2U and I4TJ]) admits a simple physical interpretation. It is the expectation value 
per unit volume £[( g ] — (Vt g \H A „ Ag g_„|O ff }/V of the flat-space Hamiltonian in the state |f2 fl ) that 
minimizes £[( g ] subject to the constraint (il g \^ a (x)\il g ) = C g ( x )i with Q (x) satisfying g a b = i] a b — 2s a b(() 
(see Eq. l5"Tj) . In other words Vl[A n , A p ,g] is the energy density in the lowest energy state subject to 
the "deformation condition" (Vt g \^ a {x)\Vlg) = ( g (x). It is precisely in this sense one should think of 
the metric g a b as determining the shape of the system. To avoid notational clutter, we will focus here 
on the case of a pure solid system and neglect the dependence on the external fields A n and A p . The 
derivation involves several steps, which we summarize below. 



First, let us evaluate the partition function in the presence of a space-time independent background 
metric g^ u of form Eq. [53] by the saddle point method. The classical solution that minimizes the 
Euclidean action and is well behaved at |a;| — > oo is given by £ a = 0. So we have: 

Z[g]=e iw ® = Exp{zVT£ (H ab (g = g,Z = Q))} (72) 

Since we are working with a diffeomorphism invariant theory, we can obtain the same result for the 
free energy in a different coordinate system. Let us use this freedom to switch from coordinates 
(x a ,g a b) to the "flat" coordinates (x a ,r] a b). The appropriate variable transformation can be 
found by noting that H ab is a scalar density. This results in x a (x) = x a — £, g (x), with the field £ g 
determined by the condition H ab (g = g, £ = 0) = H ab (g = 77, £ = £ 3 ), which explicitly reads 

W-ffl + SM«). (73) 



\dx a dib dii dx 



Eq. 1731 defines £^ up to rigid rotations and translations. For constant g a b the solution has the form 
£g (i) = K£x b +c a where the elements of K£ and c a are constant^. Equivalently the inverse change 
of variables reads x a (x) = x a + ( g L (x), with (g(x) = £, g (x) + 0(^g), and one has g a b = Vab — 2s a b((), 
with the strain s a &(C) given in Eq. 1511 

In summary, as a consequence of general coordinate invariance one has: 

Z[g] = e lW ^ = Exp [iVTC a (H ab {g = 7?,£ = Q)) } , (74) 

with a time- independent field configuration C g (x) determined by Eqs. [51] or alternatively 1731 

Next we note that the exponent on the RHS of Eq. [74] is the flat-space action evaluated at the 
field £ = ( g . Moreover, to leading order in the loop expansion (and low-energy expansion) the 
action coincides with the quantum effective action r[£ g ] = ^offfCg]. But the quantum effective 
action T[£ g ] admits an energy interpretation [36 , 37, 38 : for time-independent field configurations 
C g (x), one has that T[( g ]/T = — (n g \H^ n ^ g g =n \£l g ) where \VL g ) is the state that minimizes the 
expectation value of the Hamiltonian under the constraint (Vt g \^ a {x)\Vt g ) = Q(x) . In equations, 
the above chain of reasoning reads 

W[g] = -VTVL[g] (75) 
W[g] = S eS [Q=r[( g ]=-T(n g \H AntA9rg=ri \n g ) (76) 
= -VT£[Q , (77) 

thus proving that Q,[g] = £[( g ], with C g related to g ab by g ab = rj ab - 2s ab {( g ). 



6 The fiat coordinates (x a ,7? a 6) play a somewhat special role: the configuration § a = corresponds to the equilibrium 
configuration in absence of external fields. In this state the body-fixed coordinate z a = x a — £° are flat (coincide with the 
laboratory coordinates). Deformations from equilibrium £ ^ induce a non-euclidean metric in the body-fixed coordinates. 

7 Note however, that one needs to avoid "large diffcomorphisms" , which are not well behaved at \x\ — > oo. Proper 
behavior at infinity can be ensured by multiplying the transformation by appropriate convergence factors that decay to 
zero at \x\ — > oo faster than any polynomial. 
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